G3: Genes, Genomes, Genetics

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

Three new species of sun orchid (Thelymitra) endemic to Aotearoa New Zealand are here described. These are T. palustris, T. scabrifolia and T. semaphora. The morphological distinctiveness of these three species has been acknowledged for decades; however, their taxonomic status has remained unresolved. Evidence from existing karyological data, recently generated DNA sequence data (LFY and ycf1) and morphological studies from historical and fresh collections are used here to support their formal description. Both, T. palustris and T. semaphora are restricted to wet habitats north of Auckland (North Island). Thelymitra scabrifolia inhabits mostly scrub, and it has a similar northern North Island distribution, but is has been found also in Manawat[a]whi / Three Kings Islands and historically in Otago (South Island). All three species are polyploids and are of conservation concern.

Brettanomyces bruxellensis is a yeast species associated with diverse fermentation environments and characterized by extensive genetic diversity, including diploid, autotriploid, and allotriploid lineages resulting from independent hybridization events. These lineages are associated with distinct ecological niches and provide a framework for studying metabolic trait evolution in complex genomes. Nitrate assimilation is a relatively uncommon trait among yeasts and has been reported in B. bruxellensis, but its distribution and evolutionary history within the species remain poorly understood. Here, we combined phenotypic characterization of 151 strains with genomic analyses of 946 whole-genome sequences to investigate nitrate assimilation. Growth assays revealed that nitrate assimilation is widespread but unevenly distributed across genetic lineages, with some populations largely retaining the trait whereas others have frequently lost it. Genomic analyses identified extensive variation affecting the nitrate assimilation gene cluster composed of YNR1, YNI1, and YNT1. Nitrate assimilation was strongly associated with both gene copy number and predicted gene functionality, with nitrate-assimilating strains generally carrying more functional copies of the cluster. Leveraging the complex genomic architecture of the species, we independently analyzed primary and acquired genomes in allotriploid lineages and uncovered contrasting evolutionary trajectories following hybridization. While nitrate assimilation genes were generally maintained in primary genomes, acquired genomes showed a higher prevalence of gene loss and predicted loss-of-function variants, revealing asymmetric dynamics between subgenomes. Altogether, our results suggest that nitrate assimilation represents an ancestral trait that has been differentially maintained across B. bruxellensis lineages through a combination of copy number variation, gene degeneration, and genome-specific evolutionary dynamics. These findings provide new insights into how genome architecture and polyploid evolution shape the maintenance and loss of metabolic traits in an industrially relevant yeast species.

High-phenolic plant genera present a major technical limitation in genomic research. Standard extraction approaches that perform reliably across diverse flora often perform poorly when applied to recalcitrant taxa, producing low DNA yield and integrity incompatible with sequencing requirements. The genus Anigozanthos (Kangaroo paws) from the family Haemodoraceae exemplifies this problem. We identified key physicochemical factors governing extraction failure in this genus and resolved them through targeted modifications to lysis chemistry and contaminant management. The resulting protocol achieved a near threefold improvement in DNA purity, substantially reducing contaminant carry over and consistently yielded high-integrity, long DNA fragments (DIN > 7) across a diverse sample set spanning cultivated and wild material across four diverse genera of Haemodoraceae. We also tested a straightforward purity assessment framework that can be implemented in any standard molecular laboratory, enabling rapid pre-submission quality assessment without the need for specialised equipment. Together these advances open a practical path to genomic characterisation of Anigozanthos that establishes a transferable model for genomic research across Australia ' s chemically complex native flora.

Perennial grains represent a promising frontier for sustainable agriculture, but breeding progress is constrained by the accessibility of genotyping and the difficulty of evaluating complex traits expressed for multiple years after establishment across heterogeneous environments. Phenomic selection may help address these challenges by using inexpensive, scalable, high-dimensional phenotypes collected early in development, although the robustness of such predictions across breeding cycles remains uncertain. Here, we compared genomic selection and phenomic selection across two breeding cycles of Thinopyrum intermedium (intermediate wheatgrass; IWG; Kernza(R)), comprising approximately 2,280 individuals from maternal half-sib families evaluated across multiple field sites and years. We constructed relationship matrices from genomic markers and early-life stage phenomic data, including seed and leaf color (HSV), CropReporter multispectral reflectance and indices, and cycle-specific hyperspectral reflectance sensors. Genomic models provided the strongest predictions on average across all field traits in both cycles. Among phenomic predictors, leaf HSV was consistently the most informative, whereas CropReporter and hyperspectral data showed lower and more trait-dependent performance and seed HSV provided little predictive value. Genomic, leaf HSV, and CropReporter models transferred across breeding cycles with little apparent loss of predictive ability relative to within-cycle validation, demonstrating that their predictive signals were not restricted to a single breeding cycle. Early-life stage leaf HSV emerged as a practical, accessible tool for germplasm thinning and early-stage prioritization in perennial breeding programs. Despite limited similarity among relationship matrices, multi-relationship-matrix models rarely improved prediction beyond the stronger constituent single-relationship-matrix model. Together, these results show that early-life stage phenomic data provide reproducible information about agronomic performance expressed years later, but that predictor complexity and data integration do not guarantee improved prediction.

Developmental regression is a severe but poorly understood complication of several neurodevelopmental disorders. In Kleefstra syndrome (KLEFS1), caused by EHMT1 haploinsufficiency, regression often emerges during adolescence or early adulthood and is frequently preceded by marked sleep disturbance. Experimental work implicating EHMT1/G9a in metabolic regulation and stress responses raises the possibility that impaired metabolic resilience contributes to this vulnerability. Here, we aimed to investigate whether altered glucose metabolism is a feature of KLEFS1 and whether it relates to clinical variability, including regression. Through [18F]FDG-PET/CT, individuals with KLEFS1 who had experienced regression (n=4) exhibited a hypometabolic brain profile, whereas one individual who had not experienced regression showed globally elevated metabolic activity. In parallel, G9a mutant flies exhibited increased baseline metabolic rate and neuronal ATP levels together with sleep fragmentation resembling the clinical phenotype. Providing flies with oxidative stress to model KLEFS1 regression further exacerbated sleep disruption and was associated with a reduction in metabolic output. Importantly, adult high sugar feeding in flies prevented oxidative stress-induced worsening of sleep and maintained metabolic stability under challenge. Together, these findings suggest that regression in KLEFS1 and associated sleep disturbances are linked to underlying metabolic vulnerability and impaired maintenance of energy homeostasis under stress.

THAP9 is a transposable element-derived gene which encodes a protein that is homologous to the active Drosophila P-element transposase (DmTNP). Both THAP9 and DmTNP possess a C-terminal domain (CTD) which is functionally uncharacterized. Sequence and structural analysis suggest that the THAP9-CTD has a novel fold which is only found in THAP9 homologs. To explore the evolutionary history and characteristics of this novel domain, exhaustive phylogenetic analysis (using MSA, structure prediction, MSTA-based clustering) was performed. THAP9-CTD homologs were more widely distributed throughout the animal kingdom in comparison to DmTNP-CTD homologs which were restricted to arthropods. Moreover, the THAP9-CTD homologs were more conserved, especially among mammals and birds and their average length increased in a class-specific manner. Comparison with the DmTNP-CTD homologs demonstrates that although their respective CTDs may have evolved independently, they both surprisingly share similar secondary structure elements consisting of three conserved helical regions made of hydrophobic residues that are predicted to make up a conserved core. The role of the respective CTDs were further investigated by creating truncation mutants lacking the CTD. Interestingly both THAP9 and DmTNP truncation mutants are still capable of DNA excision and integration suggesting that their respective CTDs are not essential for DNA transposition. Moreover, CTD truncation favours DNA integration in THAP9: this suggests that CTD acquisition during evolution may have led to THAP9 domestication as observed in other transposable element-derived genes like Rag1 and piggybac, which have similar terminal regulatory domains.

Myxococcus xanthus bacteria form aggregates when starved on solid surfaces and some cells differentiate into spores. Studies of mutants in monoculture have advanced knowledge of this multi-cellular developmental process, but our understanding of the genetic determinants is incomplete. To assess gene function genomewide, we generated a pool of barcoded transposon insertion mutants, subjected it to starvation, and separated developmental samples into non-aggregated cells, aggregated cells, and spores. We also subjected our pool to chemically-induced unicellular sporulation. Evaluation of changes in the abundance of mutants in samples allowed identification of 200 genes in which insertions reproducibly caused distinct patterns of depletion and/or accumulation over time. Many of these genes have well-established roles in development, validating our approach, while many others have not previously been associated with development. Genes involved in type IV pili (T4P)-dependent motility were more important than gliding motility genes for aggregation and sporulation in the mixed population. Although exopolysaccharide (EPS) synthesis genes are required for aggregation in monoculture, most were dispensable for aggregation in our pool, consistent with EPS sharing between cells, yet these genes were required cell-autonomously for efficient sporulation. Genes for positive regulators of EPS synthesis were important for aggregation as well as sporulation, suggesting functions beyond EPS production. Insertions in several novel genes impaired both starvation- and chemically-induced sporulation. Many genes increased the efficiency of starvation-induced sporulation. Some of these mutants, which we call "developmental winners", are novel cheaters. Our results demonstrate the power of using the newly-created mutant library to elucidate M. xanthus biology.

Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.

Background: Reproducible taxonomic collapsing and geological-timescale annotation of time-calibrated phylogenetic trees in R often require coordination among several packages and repeated code for label parsing, clade validation, plotting, and export. Workflow-managed analyses additionally benefit from non-interactive configuration, predictable diagnostics, and machine-readable exit status. Results: We present Rclade, an R package that consolidates the multi-package coordination required for taxonomic collapsing into a streamlined, single-function interface. Rclade provides (1) custom ggproto objects (GeomPolygonStraight/GeomSegmentStraight) that bypass coord_munch() interpolation to achieve straight-edge rendering of collapsed triangles in circular layouts; (2) automatic detection and parsing of four taxonomic-label formats (GTDB, Silva, NCBI, embedded) plus user-supplied custom regex, with explicit input-validation contracts and parsing-accuracy evaluation on real and derived test sets; and (3) workflow embeddability through YAML configuration, library-mode APIs, and standard Unix exit codes. Benchmarks on synthetic and real datasets (200-10,000 synthetic tips and real reference trees up to 10,122 tips; 5 replicates at every scale under a unified fully rendered measurement protocol) show that the full-pipeline overhead is modest for interactive use (median {approx}0.87 s in-session rendering and {approx}8.4 s process-level wall-clock at 10,000 tips). Conclusions: Rclade is a convenience layer over the ggtree/deeptime ecosystem that reduces boilerplate while adding targeted technical improvements for circular-layout rendering and format heterogeneity management.

Nuclear pore complexes (NPCs) are nuclear envelope (NE)-embedded protein assemblies that mediate nucleocytoplasmic exchange and interact with the genome, including binding of an NPC component Nup93 to Polycomb chromatin domains. Here, we investigated the in vivo relevance of this relationship in Drosophila, which unusually contains two distinct paralogs of Nup93. Interestingly, we identified a Nup93-2-specific tumorigenic phenotype in larval wings, where depletion of Nup93-2, but not Nup93-1, led to tumor-like overgrowth, reminiscent of Polycomb mutations. Consistently, our transcriptomic analysis revealed a wide-spread loss of gene silencing in Nup93-2-depleted wings, particularly in a Nup93-bound Polycomb domain spanning genes for activators of JAK/STAT signaling. Nup93 paralogs were not found to differ in their effect on NPC biogenesis but strikingly, showed differences in subnuclear localization patterns. While Nup93-1 co-localized exclusively with fully assembled NPCs, Nup93-2 exhibited only partial co-localization and was found at additional NE locations in a tissue-specific manner. Together, our results identify an in vivo silencing role of a Nup93 paralog and suggest that Nup93-2 may form a unique NE-associated complex that targets a subset of Polycomb domains containing growth-promoting genes.

Bacterial genome-wide association studies (GWAS) offer a powerful approach to identify the genetic basis of a trait measured in a set of sequenced isolates. As the number of sequenced isolates has grown, the limiting factor for GWAS has become phenotyping enough isolates to achieve statistical power. To overcome the need for large-scale phenotyping, we developed Bayesian Adaptive Sequential Sampling GWAS (BASS-GWAS), which couples Bayesian adaptive experimental design with a sparse regression model to select maximally informative isolates for phenotypic testing. BASS-GWAS efficiently recovered causal loci for three antimicrobial resistance traits in Neisseria gonorrhoeae, requiring many fewer phenotyped isolates than random sampling. We applied BASS-GWAS to discover variants enabling gyrBD429N-dependent cross-resistance to the novel topoisomerase inhibitors zoliflodacin and gepotidacin. After phenotyping fewer than 30 isolates, we identified and then validated both parCD86N and a gyrA-parE-based pathway as enabling cross-resistance. BASS-GWAS provides a practical and statistically principled solution for efficient bacterial GWAS.

Pre-mRNA splicing is essential for gene expression, yet how disruption of core spliceosomal factors produces tissue- and developmental stage-specific phenotypes remains poorly understood. Here, we investigated the in vivo function of the conserved spliceosomal kinase PRPF-4 in C. elegans using endogenous reporter analysis, conditional protein depletion, and transcriptome-wide analysis of alternative splicing and gene expression. We found that PRPF-4 is broadly expressed throughout development and is continuously required for postembryonic development, with distinct requirements in the pharynx, nervous system, and germline. Acute PRPF-4 depletion rapidly disrupts alternative splicing across thousands of transcripts, with exon skipping representing the predominant class of affected events. In addition, PRPF-4 depletion results in a robust transcriptome shift with induction of components of the spliceosome and repression of ciliary and ion transport-related transcripts. These findings establish PRPF-4 as a central regulator of RNA metabolism and demonstrate the far-reaching effects on gene expression caused by loss of core spliceosomal components.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

MicroRNAs (miRNAs) are key post-transcriptional regulators of development, stress response, and secondary cell wall formation in woody plants, yet annotations for Eucalyptus grandis, the world's most widely planted hardwood, remain fragmented across studies using incompatible discovery pipelines and filtering criteria. Here we present the Eucalyptus MicroRNA Archive (EMA), a curated, locus-resolved database integrating three independent small RNA sequencing datasets spanning vegetative tissue, somatic embryogenesis, and mechanically induced tension wood formation. Applying annotation criteria aligned with current plant miRNA standards, EMA catalogs 99 curated miRNAs (31 previously described, 68 novel) organized into 34 family-level groupings under a three-tier confidence system, known-reference-supported, multi-study replicated, or single-study, that preserves study-of-origin and sample-level evidence for every entry. Cross-study comparison showed that only 9 of 99 entries (9.1%) were independently supported by all three datasets, supporting an evidence-tiered rather than binary annotation scheme. Target prediction against the E. grandis transcriptome yielded 1,773 miRNA-target interactions spanning 764 loci, integrated into a combined miRNA-target and protein-protein interaction network. This network resolved into functionally coherent, mutually isolated clusters, including an miR482-associated NBS-LRR/TIR disease-resistance hub with a substantial translational-repression component, alongside modules enriched for ribosome biogenesis and translation, DNA replication, and nitrogen and carbohydrate metabolism. EMA is publicly accessible through an interactive web dashboard, with all curated data, source code, and analysis scripts openly available, providing a reproducible, extensible framework for E. grandis miRNA research and a template for similarly structured resources in other non-model woody species.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

The exponential increase in the number of genotyped animals, combined with the availability of high-density SNP chips has introduced computational challenges for routine genomic evaluations, particularly during the construction of the genomic relationship matrix. Although higher-density SNP panels can facilitate the identification of causal mutations, their use substantially increases computational requirements without a proportional gain in genomic prediction performance. To optimize computational efficiency while maintaining accuracy of genomic predictions, this study compared five SNP selection strategies (i.e., random sampling, random sampling with inclusion of informative SNPs, linkage disequilibrium (LD)-based pruning, a Shannon entropy-based machine learning approach, and [[EQUATION]]-based prioritization) to develop reduced-density panels for Nellore cattle. Using high-density (HD) genotype data comprising 437,650 SNPs from 304,782 animals (after quality control) as reference, three reduced-density panels (25K, 45K, and 65K SNPs) panels were tested across five traits (i.e., Age at first calving, Stayability, Weaning weight, Yearling weight, Muscling) with diverse genetic architectures. Genomic estimated breeding values (GEBVs) derived from these reduced panels were compared to those obtained from the HD reference panel using Pearsons correlations, under both genomic best linear unbiased prediction (GBLUP) and single-step GBLUP (ssGBLUP) methods. In the GBLUP model, prediction accuracy generally improved with increased marker density. Random selection with and without the informative SNPs consistently yielded the highest accuracies, whereas the [[EQUATION]]-based approach showed the lowest agreement with the HD reference across all densities. In contrast, ssGBLUP demonstrated strong robustness to marker reduction, producing uniformly high correlations {approx}1.00) across all SNP densities and selection strategies. These findings indicate that optimized low-density SNP panels maintain prediction accuracy comparable to HD panels, offering a cost-effective tool for large-scale genomic evaluations.

Thermodynamic models quantitatively describe interactions between transcription machinery and bacterial promoters. Contrary to conventional understanding, model analysis by Parisutham et al. (2025) attributes transcriptional inhibition by repressors to overstabilization of the RNA polymerase-promoter complex rather than prevention of its formation. Moreover, it suggests an inverse scaling relationship between basal promoter strength and transcriptional fold change, applicable to both repressor- and activator-mediated regulation. To reevaluate findings from this study, we systematically analyze empirical data and compare its framework with conventional thermodynamic models. In contrast to the inverse scaling relationship, data across multiple sources exhibit a peaked tradeoff between basal promoter strength and fold change, underscoring the importance of broad data coverage in revealing the full pattern required for reliable model inference. Furthermore, we identify the model assumption responsible for the apparent inverse scaling and misinterpretation of regulatory mechanisms. Relaxing this assumption enables the model to capture the peaked tradeoff and yield inferences consistent with established mechanisms of transcriptional repression and activation. We further derive a mathematical solution that connects basal expression to fold change for both repressor- and activator-regulated promoters. Our results underscore the importance of broad data coverage to avoid a blind-men-and-elephant interpretation and establish basal promoter strength as a key design parameter governing transcriptional regulation.

Elucidating the biological and molecular mechanisms that govern male fertility and mating behavior in mosquitoes is critical for optimizing genetic and sterile insect technique-based vector control strategies. Here, we examined age-related changes in male reproductive capacity in Anopheles culicifacies, using female egg output as an indirect indicator of male fertility. Our results demonstrated that male reproductive age follows a non-linear pattern of fertility. Morphometric analysis from emergence to day 13 post-eclosion revealed a strong correlation between seminal vesicle capacity and female fecundity, suggesting that age-dependent gonadal development directly influences reproductive potential. At the molecular level, we identified AcLITAF6 as a key regulator of male reproductive homeostasis. RNAi-mediated knockdown of AcLITAF6 impaired apoptosis-associated and phagocytic clearance, reduced sperm viability, and decreased female productive outcomes. Conclusively, we reveal a previously unrecognized role of LITAF in sperm quality control and male reproductive fitness, highlighting AcLITAF6 as a potential target for mosquito population suppression strategies.

Bovine viral diarrhea virus (BVDV) remains an economically important pathogen of cattle despite widespread vaccination. A homozygous CD46-edited Gir heifer (Ginger) was previously shown to have significantly reduced susceptibility to BVDV. The edited allele contains an in-frame six amino acid substitution within the virus-binding domain of the BVDV entry receptor CD46, replacing residues G82QVLAL with A82LPTFS. Here, we investigated whether reduced BVDV susceptibility is maintained when the edited allele is inherited in the heterozygous state. Ginger was artificially inseminated with semen from an unedited Gir bull and produced a healthy heterozygous CD46-edited bull calf (Giraldo). Whole-genome sequencing confirmed the inheritance and structural integrity of Giraldo's edited allele. Compared with Ginger, Giraldo exhibited similarly reduced ex vivo BVDV susceptibility across primary fibroblasts, lymphocytes, and monocytes, despite inheriting a wild-type CD46 allele from the sire. Allele-specific CD46 RNA expression analysis demonstrated expression of both the edited and wild-type CD46 alleles. Thus, the reduced-susceptibility phenotype was not attributable to transcriptional silencing of the wild-type allele. Lentiviral complementation studies in CD46-knockout Madin-Darby bovine kidney (MDBK) cells further demonstrated that this wild-type CD46 allele was competent to support BVDV infection when expressed independently. Together, these findings indicate that the CD46 A82LPTFS allele can confer reduced BVDV susceptibility in the heterozygous state despite expression of a functional wild-type CD46 allele. This result suggests the potential to more rapidly disseminate reduced BVDV susceptibility through conventional breeding using homozygous CD46-edited sires.